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I Polarization observations of the cosmic microwave background with the Cosmic Back- 

ground Imager from September 2002 to May 2004 provide a significant detection of 
the E-mode polarization and reveal an angular power spectrum of polarized emission 
showing peaks and valleys that are shifted in phase by half a cycle relative to those of 
the total intensity spectrum. This key agreement between the phase of the observed 
polarization spectrum and that predicted on the basis of the total intensity spectrum 
provides support for the standard model of cosmology, in which dark matter and dark 
energy are the dominant constituents, the geometry is close to flat, and primordial 
density fluctuations are predominantly adiabatic with a matter power spectrum com- 
mensurate with inflationary cosmological models. 

In recent years a wide variety of observations have provided support for a standard model of cosmology 
and cosmic structure formation. In this model {1 2 3), the mass-energy density of the universe is domi- 
nated by cold dark matter and dark energy, possibly in the form of Einstein's cosmological constant, and 
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conventional baryonic matter is only a minor component. As the universe expanded from its hot, dense 
origins, all the structures seen in the universe today formed under the action of gravity on initial nearly 
scale-invariant adiabatic gaussian density fluctuations. Observations of the anisotropy in the intensity of 
the cosmic microwave background radiation (CMB) have provided much of the evidence for this model 
and estimates of the values of the fundamental parameters, including the spatial curvature and the den- 
sities of dark energy, cold dark matter, and ordinary matter. The two major ingredients of the standard 
model, dark matter and dark energy, are far from understood, and their existence presents the most seri- 
ous challenge to physics since the quantum and relativistic revolutions of a century ago. It is therefore 
essential to extract as much information as possible from the observations in order to test all aspects of 
the standard model and to look for possible anomalies that might provide insights into the nature of these 
two dark components of the cosmos. 

In the 1980s successively more stringent limits were placed on the observed level of temperature 
anisotropy in the CMB {4 5), providing convincing evidence that the dominant matter constituent in 
the universe is non-baryonic. These searches culminated in the detection of anisotropics by the CORE 
satellite @, a confirmation of one of the major predictions of theoretical cosmology (ITI 151 191 17^ 1771 
17^ . Rapid advances in experimental techniques have since delineated the prominent features in the 
angular power spectrum U3.J4 15 16.17 18). The spectrum of fluctuations on large angular scales [low 
multipole numbers, / < 500 ( 1791 11 has been measured with high precision by the Wilkinson Microwave 
Anisotropy Probe (WMAP) STUi . whereas precise and sensitive CMB observations from the ground and 
from balloon-based platforms have extended the spectrum to angular scales as small as a few minutes of 
arc {I w 3500) QL,22ii23^^^5). 

The early universe was opaque to electromagnetic radiation, but as it expanded and cooled the hot 
electron-baryon plasma combined into neutral hydrogen and helium and the universe became transparent. 
The microwave background photons that we detect today have passed freely through the universe since 
they were last scattered by electrons in the ionized plasma. They thus provide a picture of the physical 
conditions at the time of last scattering when the universe was about 400,000 years old. The angular 
power spectrum of the CMB reveals the initial fluctuation spectrum modulated by the effects of acoustic 
waves in the plasma i26i . and it gives quantitative information about the physical conditions in the 
plasma. 

The polarization of the CMB provides an independent way to test the standard model (I27II2^|B91 I91 
1^ . Anisotropic Thomson scattering of photons at the time of last scattering gives rise to weak linear 
polarization of the CMB. Measurement of the CMB polarization not only provides an additional way to 
measure the parameters of the model, it can also verify the correctness of many of the basic assumptions 
on which the model is founded. It is for this reason that many experiments are being designed to measure 
the CMB polarization power spectra, despite the difficulty of the observations. 

Angular power spectra give the variance C/ (usually expressed in terms of CMB temperature, and 
with units of //K^) as a function of multipole number I (l26l l. The intensity of polarized radiation can be 
expressed by the four Stokes parameters /, Q, U, and V i31\\32l . Total intensity is represented by /, 
linear polarization by Q and U, and circular polarization by V. Thomson scattering does not generate 
circular polarization, so we ignore V. From the other parameters we can generate three power spectra 
TT, EE, and BB, where T is the total intensity (Stokes /) and E and B are the curl-free and curl-like 
components of the linear polarization field (Stokes Q and U) i33l34l35l and also three cross-spectra TE, 
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TB, and EB. Because of the parity properties of the T, E, and B signals, the only non-zero spectra 
should be TT, EE, BB, and TE. In the standard model £'-modes are generated from the primary scalar 
density fluctuations, whereas S-modes are generated only by gravitational-wave tensor fluctuations and 
secondary processes; the predicted BB power spectrum is undetectable with cun^ent sensitivity. A high 
level of BB would require modifications to the standard model, or it could indicate that the observations 
are contaminated by radiation from foreground sources, because these are epected to produce both E and 
B modes in equal measure. The polarization spectrum of the CMB is more difficult to study than the total 
intensity spectrum because the fractional polarization of the CMB radiation is no more than 10%. After 
a number of experiments that placed upper limits on CMB polarization, EE power has been detected by 
the DAS I [Degree Angular Scale interferometer] experiment (with 6.3 cr significance) i36l37}\38l and the 
CAPMAP [Cosmic Anisotropy Polarization Mapper] experiment {2.3 a) OU- TE cross-spectral power 
has been detected by DASI and by WMAP {40). 

The TT spectrum arises from density and temperature fluctuations in the plasma, but polarized radi- 
ation, which is caused by the local quadrupole at the time of last scattering, is sensitive to the velocity of 
the plasma. Because velocity and density are out of phase in an acoustic wave, the maxima in the EE 
spectrum are out of phase with those in the TT spectra. This phase shift between the spectra is a key fea- 
ture of the standard model. It has been seen at large angular scales in the TE spectrum by WMAP {47^, 
but it has not yet been verified directly through EE or at the small angular scales corresponding to 
clusters of galaxies. 

We report here observations made with the Cosmic Background Imager (CBI) that have sufficient 
sensitivity and resolution to detect and measure the second, third, and fourth peaks in the EE spectrum 
i42i . determine the TT-to-EE phase shift, and thus further test the standard model. 

The Cosmic Background Imager 

The CBI (Fig. [0 has been making observations of the CMB from a site at 5000 m elevation on the 
Chajnantor plateau in the Chilean Andes since late 1999. It is a 13-element radio interferometer receiving 
radiation in 10 I-GHz frequency channels covering 26 to 36 GHz i43\\44l . The individual antennas are 
0.9 m in diameter, and the possible baselines range in length from 1.0 to 5.5 m. An interferometer 
baseline of length d is sensitive to multipoles / around 2TTd/X where A is the observing wavelength. The 
CBI can thus measure the spectrum from / 300 to / ~ 3500. The antennas are mounted on a platform 
with azimuth and elevation axes that allow all the antennas to track a point on the sky. The platform can 
also be rotated about the line of sight; this allows full sampling of all possible baseline orientations and 
facilitates calibration of the instrumental polarization effects. 

Each antenna is sensitive to a single sense of circular polarization, right (R) or left (L). Co-polar 
basehnes, RR and LL, are sensitive to Stokes I ±V ^ I (assuming circular polarization is negligible), 
whereas cross-polar baselines, RL and LR, are sensitive to linear polarization. Stokes Q ± ill (45). It 
was thus straightforward to adapt the CBI to measure linear polarization by changing the sense of some 
of the antennae to maximize the number of cross-polar baselines [Supporting Online Material (SOM) 
Text]. 

The observations reported here were carried out between 22 September 2002 and 7 May 2004 using 
7 antennae with polarizers set to left circular polarization (L), and 6 antennae with polarizers set to right 
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circular polarization (R). To avoid contamination by the sun and moon, observations were made only at 
night and at angles of greater than 60° from the moon (supporting online text). 

The size of the CBI antennae sets the / resolution of observations made in a single pointing to Al w 
300. To obtain the higher resolution in / necessaiy to resolve the expected structure in the EE spectrum, 
we had to image a larger area by making a mosaic of overlapping pointings. From 2002 to 2004, we 
observed a grid of pointings in four regions near the celestial equator that were separated by about 
6 hours in right ascension and identified as the 02^, 08^, U^, and 20"^ fields (Fig. 01. These fields 
were centered on those we observed in 2000 and 2001 in order to measure the TT power spectrum 
( l?Hll?7ll?^l?Pll5Dll5710?t . The sepai^ation of the pointings was 45 arcmin, twice that of the earlier 
observations, leading to modulation of the sensitivity across the field. For three of the fields we used 
36 different pointing positions giving fields « 5° square, but for the 20"^ field we divided the available 
integration time between six pointings in a row; these deeper observations should be more sensitive to 
any potential systematics. 

The largest source of diffuse foreground contamination over the 26- to 36-GHz band is synchrotron 
radiation from the Galaxy. We chose the four CBI fields, which were constrained to be separated by 
about 6 hours in right ascension, so as to minimize this contamination (Fig.l^Jl. The CBI 02^, 08^, and 
20"^ fields, like the DASI polarization field, are in regions of low synchrotron emission (Fig. but the 
14'' field is near the North Polar Spur, and the WMAP observations suggest that this has a higher level 
of synchrotron foreground. The wide separation of the CBI fields provides some control on foreground 
contamination, because foreground emission is unlikely to be correlated over such large distances and if 
it were a problem we would expect to see differences between the spectra of the different fields. 

The largest systematic instrumental effect we have to eliminate is ground spillover. Although the 
ground radiation is unpolarized it enters the CBI feeds after reflection off the inner surface of the shield 
cans that reduce cross-talk between receivers and therefore gives rise to a highly polarized contam- 
inating signal. This signal is particularly strong on the shortest baselines at the lowest frequencies. To 
eliminate the ground radiation, we observed sets of six fields separated by 3 min in right ascension at the 
same azimuths and elevations, spending 3 min on each field, so that if the ground emission is constant it 
should make equal contributions to all six fields. When estimating the power spectrum we made use only 
of the differences between the fields, ignoring the contribution that is common to all six. This strategy 
requires that the ground be stable over the total scan duration of 18 minutes, but the penalty is only ^6/5 
in flux density sensitivity, equivalent to a factor of 1 .2 in observing time. 

The data were edited to remove data corrupted by instrumental or atmospheric problems. Amplitude 
calibration was based on Jupiter (14) and polarization position-angle calibration was based on Tau A, for 
which we measured a polarization position angle {E vector) of —27.6° by comparison of CBI and Very 
Large Array ) observations of 3C 273 and 3C 279. Instrumental polarization leakage was measured 
on Tau A and found to be negligible. The noise was estimated from the scatter of the measurements in 
each 18-minute scan (SOM Text). 

After the editing and calibration of the data we made images of the four fields and of all the calibration 
sources in order to check for possible anomalies {48). As an example, the / image of the field (Fig.ill, 
made without any subtraction of ground spillover or foreground sources, shows significant power above 
the level of the noise. This is due to both CMB emission and ground spillover. In the Q and U images 
(Fig. the signal from the CMB is too weak to identify, and these images are dominated by the regular 
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pattern due to ground spillover. The level of the ground spillover in Q and U indicates that there is 
some ground contamination in the / images as well, although it is somewhat weaker than the CMB 
signal. When we estimate power spectra, the ground spillover is removed from the data by projecting 
out (53 21 ) the common mode in the six matched pointings, so the visibility data set from which the 
images of Fig. |4] were made is the data which we use in the CMB spectrum determination. But in 
order to check our procedures we have also made images from the differences of visibilities measured 
in pointings separated by 9 min in right ascension; these images should be free of ground contamination 
(Fig- El- The total intensity / image shows power well in excess of the noise level, whereas the Q and U 
images show only noise, the sensitivities per resolution element being too low to reveal the polarization 
of the CMB. These images also show that leakage of total intensity into the polarization data is small 
compared to the noise. 

Power Spectrum Estimation 

To estimate power spectra from the interferometer visibility measurements we used maximum-likelihood 
procedures similar to those adopted for earlier experiments (15413 7II55I) . To process the CBI data, we have 
extended the gridding-based procedure used in our earlier work to deal with mosaicked polarization 
observations. A given correlator output sample, or visibility, can be one of the four polarization products 
RR, RL, LR, or LL. These can be related to the fundamental CMB polarization modes T (temperature), 
E, and B (polarization) {37). The covariances between the measurements depend on the six CMB 
covariances TT, EE, BB, TE, TB, and EB. Because the CBI measures circular polarization products, 
which are orientation independent (depending only on the handedness of the wave polarization), the CBI 
(or any interferometer using circularly polarized receptors) is sensitive to the E and B modes directly. 
This simplifies the power spectrum analysis (SOM text). 

The principal foreground contamination in total intensity for the CBI is that due to extragalactic 
radio sources WA. For the total intensity spectrum, in which discrete sources have a substantial impact, 
our approach is similar to that used in earlier CBI analyses ( 47.. .48... 24.. 49) . with minor modifications. 
Some 3727 NRAO Verly Large Array (VLA) Sky Survey (NVSS) ^ sources with SiagHz > 3.4 mJy 
were projected out of the data. In previous work we used separate covariance matrices with different 
projection factors for sources that were detected at 32 GHz on the OVRO 40-meter telescope and for 
those that were not. Because in the end there was no gain from this approach, in the present analysis we 
combined all sources into a single covariance matrix with a single projection factor. For this analysis we 
assume a uniform variance of 1 Jy^ for each source, rather than adjusting the variance for each source. 
We find that this yields matrices that are numerically more stable under the action of our procedure of 
completely projecting the source modes out of the data. After a number of tests, we adopted a value of 
gsrc = 100 for the pre-factor (equivalent to setting the variance on each source flux density to 100 Jy^). 

However, because non-thermal extragalactic radio sources are weakly polarized (P < 10%) and 
furthermore only a small fraction of them have P > 2%, only a few of the sources that we projected out 
in the total intensity spectrum can affect the polarization spectrum. When estimating EE and BB, we 
therefore projected out only a subset of the NVSS sources with 5i.4 > 3.4 mJy. In total 556 of these 
potentially troublesome point sources need to be considered. These include (i) NVSS sources with > 3a 
detections of polarized flux density at 1.4 GHz and (ii) sources detected by the 30-GHz OVRO survey of 
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the 2000-to-2001 CBI total intensity fields (W). The projection of 556 sources out of the CBI data has 
only a small effect on the EE power spectrum: in all bins the effect is <^ Icr. In the first two bins, where 
the polarization detection is strongest, the effect is less than 3 /iK^ for each bin. Both with and without 
projection, the BB power spectrum is consistent with zero, and the EE spectrum changes very little. 
Therefore, the sources that we have identified as potentially troublesome (with the criteria described 
above) have a negligible effect. If we have failed to identify some sources (highly polarized sources 
just below the NVSS detection limit, for instance), they should add a characteristic /(/ + 1)C/ oc P 
contribution to both EE and BB power spectra, and show up more strongly in the lower frequency 
channels. No such signature is evident in the CBI data. 

We have also studied the effects of point sources using Monte Carlo simulations. To do so we 
used the NVSS source statistics to characterize the fractional linear polarization of sources, finding a 
mean 1.4 GHz polarization of 2.7%. Most sources had polarizations less than this; 4% had polarizations 
greater than 10%, and 1% had polarizations greater than 15%. Because the fractional polarization of 
CMB anisotropics is ~ 10%, the discrete source foreground will be relatively weaker in polarization 
than in total intensity. In the EE power spectrum analysis of the simulated data, we find that the first 
two bins change by less than 4 /iK^ when the source projection is turned on, similar to what is seen in 
the real CBI data. 

The common-mode signal from the ground was removed by constructing a scan covariance matrix 
assuming a unity correlation between identical visibilities coming from the same scan (e.g., for the six 
visibilities taken in the consecutive 3 minute integrations that constitute a scan), then passed through the 
gridding operation. The modes defined by this scan covariance matrix were projected out of the data ( 53j 
by applying a large prefactor to this matrix in the likelihood maximization procedure (essentially setting 
the variance of these modes to be infinite), in the same manner as for the point sources. We used simulated 
data to determine the best value of the prefactor, and found that too small a value did not completely 
eliminate the ground spillover, whereas too large a value caused numerical problems. Because there are 
a large number of sources in the list used for the TT projection, there is an interaction between the source 
and scan projection mati^ices gsrc C'^'''^ and qscan C^^^'^ when the pre-factors Qsrc and gscan become large. 
We explored a range of values for these, and found that for the value of ggrc = 100 a value (/scan = 100 
was in the center of the range for which the TT band powers were stable (there was no significant change 
in TT amplitude from gscan = 10 to Qscan = 100). Similar tests on the real data also showed that the 
ground signal was eliminated whereas the band powers remained stable. 

Polarization Power Spectra 

The CBI measurements for all 10 frequency channels and all four fields observed from 2002 to 2004 
have been combined in the maximum likelihood procedure to estimate the TT, EE, TE and BB power 
spectra (Tableland Fig.|6ll. The scan means (for ground contamination) and point sources have been 
projected out as described in the previous section. We divided the I range into seven bands, with most of 
the bins having width Al = 150. Adjacent bands are anticorrelated at the 10 to 20% level (fig.|9ll. Finer 
binning is possible, but this gives larger band-to-band correlations and is less satisfactory for presenta- 
tion. For the quantitative analysis below, we have used bins of width Al « 75 and taken into account 
the band-to-band correlations. The results from both binnings are consistent. We have also calculated 
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window functions tliat can be used to calculate the expected band powers in our bands from a theoretical 
model spectrum (fig.fTOb. 

We compare our results with a fiducial model spectrum (Fig. |6ll, for which we have chosen the 
theoretical model C57j with a power law for the primordial spectral index which best fits the first-yeai" 
WMAP, 2000 CBI, and ACBAR [Arcminute Cosmology Bolometer Array Receiver] CMB total-intensity 
data [the "WMAPext" data set (T)]. Our results are consistent with the predictions of this model. We 
have checked this by calculating for a comparison of our measured band powers and the band powers 
predicted by the model, with the CBI window functions and the full band-to-band covariance (estimated 
from the Fisher matrix). The values of [for 7 degrees of freedom (df)], with the probabilities of 
obtaining larger values under the null hypothesis in parentheses, are: for TT, 7.98 (probability = 0.33); 
for EE, 3.77 (probabihty = 0.80); for BB, 4.33 (probability = 0.74); and for TE, 5.80 (probabihty = 
0.56). 

The TT spectrum shows the same features that we saw in the CBI 2000-to-2001 observations 
the most prominent being the drop in power between the third and fourth acoustic peaks. The TT 
spectrum from 2002-2004 is slightly higher than the fiducial model, but the difference is not significant. 
Both the EE and TE spectra are consistent with the predictions of the fiducial model. The EE spectrum 
shows detection of power at / < 800, whereas the TE spectrum is not sufficiently sensitive to show a 
positive detection. No power is detected in the BB spectrum, as expected on the standard model. The 
95% confidence upper limit on BB power (assuming flat band power in a single I bin) is 7.1 /xK^ 
with source projection or 2.7 /xK^ without source projection. Because ground radiation and foreground 
sources are expected to contribute equally to EE and BB, this low limit on a possible BB component 
at multipoles I < 1000 demonstrates that there is no significant ground or point source contamination at 
these multipoles and gives confidence in the reliability of the EE spectrum. 

At present, the addition of polarization data to CMB TT data has little effect on the values and 
precision of cosmological parameter estimates, because of the weakness of the polarized signal relative 
to the total intensity signal. Rather, the strength of the measurement of EE lies in its ability to test a 
different aspect of the theory. It is nonetheless interesting to explore the effect of the new polarization 
results on cosmological parameter estimation and to check consistency. As sensitivities improve, future 
polarization data will have a bigger impact on parameters (JS'). In this initial investigation including CMB 
EE polarization data, we explore a limited set of cosmological parameters that has been successful in 
describing all aspects of CMB data. The model has its basis in the simplest inflationary paradigm, 
characterized by the following basic set of six parameters: loi^ = Qi^h?, the physical density of baryons; 
iOc = r^c/i^, the physical density of cold dark matter; = 100^7^, parameterizing the angular scale 
associated with sound crossing at decoupling, which defines the overall position of the peak-dip pattern; 
Us, the spectral index of the scalar perturbations; In A^, the logarithm of the overall scalar perturbation 
amplitude; and t^, the Thomson scattering depth to decoupling. We do not consider any gravitational- 
wave induced components because they are not expected to be detectable by the CBI. 

The strongest prior we impose is that we only consider flat models (l^tot = 1)> as expected in most 
inflation models. We note that the parameters we derive can change significantly when this prior is 
relaxed We also impose a weak-/i prior comprising limitations on the parameter h (0.45 < h < 

0.90, where h = i7o/(100 km s^^ Mpc^^) and Hq is the Hubble constant) and on the age of the universe 
(to > 10 Gyr). Within the context of flat models the weak-/i prior influences the results very little. We 
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do note, however that extreme models with high Thomson depth are excluded by this prior. 

In our analysis we consider three combinations of data: (i) WMAPl only TT and TE results from the 
first year WMAP data (40), using the likelihood procedure described in (60); (ii) CBI pol plus WMAPl 
obtained by adding 14 band powers in each of TT, TE, and EE obtained from an analysis of the 2002- 
to-2004 CBI data with a bin width A/ w 75; and (iii) CBI pol plus CBIext plus WMAPl, consisting of 
all of the CBI polarization data and the addition of high-/ bands from our combined mosaic and deep 
field TT results (24 ), covering the range / = 600 to / = 1960 [bands 5 to 14 of table 1 of (Ell ■ This third 
combination extends the data well into the region known as the damping tail, where power is suppressed 
by photon diffusion and the finite thickness of the last scattering surface. 

We use a modified version of the Markov chain Monte Carlo (MCMC) package COSMOMC (I67l62t 
to evaluate probability distributions of the various parameters with respect to the CMB data. We have 
extended our earlier procedures (I24ll59t to include polarization spectra and the cross correlation between 
TT and EE spectra (63). In addition to estimating the six cosmological parameters defined above, we 
determined the distributions of six other derived parameters from the same Markov chains: I^a the energy 
density in a cosmological constant in units of critical density, the total age of the universe in Gyear; the 
total energy density of matter, fi^; the present-day RMS mass fluctuation on 8/i~^Mpc scales, cxg; 
the redshift of reionization, Zre (related to and il^); and lastly the Hubble parameter Hq in units of 
km s^^ Mpc^^. is a derived quantity determined from 6. The amplitude parameter cjg is related to 
In ^4^ and has more relevance for comparison with large scale structure data. As expected the inclusion 
of our polarization results does not have a large impact for this limited parameter set (Table|2li. However 
when including the CBIext TT band powers we obtain significant reduction in the uncertainties, in 
agreement with (.24 J . 

Significance of detection. Our standard CMB power spectrum analysis d6?l involves the use of a 
fiducial Q shape against which the band powers are evaluated. The gridding procedure breaks the power 
spectrum into top-hat bands in and thus the (multiplicative) band powers qb effectively break the 
spectrum into piecewise continous bands that follow the shape Cbi = Qb Ci for / within each band B. 
The most conservative choice for a shape is the flat spectrum Ci = 2it/P, but one can use a matched 
shape derived from an actual CMB power spectrum and thus optimally check for deviations from that 
model. This also allows the use of wider / bands. If we use the fiducial model fitted to the WMAPext 
dataset (Fig.l^Jl as our shape, and only project point sources from the TT sector, we find for the CBI data 
in a single I band a maximum likelihood value band power for EE of qs = 1-22 it 0.21 (68%) with 
respect to the WMAP-normalized spectrum, with a value for the log-likelihood with respect to zero of 
39.8 (equivalent to an 8.9 a detection, where a = \J2 A log £). This can be compared with the detection 
of 6.3 a reported for the DASI 3-year results Although there is no indication that the polarization 
of the foreground point sources is affecting our data, we can also adopt a conservative approach and 
project out the subset of the brightest sources, as described above. In that case, the best-fit band power 
qB is 1.18 lb 0.24 (68%) with log-likelihood with respect to zero of 24.3 (equivalent to T.Oct). This 
reduction in significance is due to the increase in uncertainties from the lost modes in this projection, 
i.e., the drop in band powers is negligible, which again suggests that point sources are not a problem in 
the EE spectrum. Although we find no evidence for point sources affecting our EE spectrum, we adopt 
this more conservative value as our estimate of the significance of our detection. 
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Phase of the acoustic oscillations. The measurement of the phase of the polarization EE spectrum 
can, in principle, provide one of the fundamental pillars of the standard model because it tests a unique 
aspect of the acoustic waves in the photon-baryon fluid. The peak positions in TT are proportional to 
irlsj, whereas for EE polarization they are proportional to TTls{j + 1/2), with some corrections from 
projection effects. To test this, we devised phenomenological models in which the phase-relationship be- 
tween TT and EE is changed. For these models, we first approximated the fiducial model EE spectrum 
as a function: 

/(/ + l)C7f ^ = f{l) + g{l) sm{kl + 0) (1) 

where / and g are smooth, non-oscillating functions (we used rational functions with quadratic numerator 
and denominator) and /c is a constant. We then varied to get a range of phase-shifted spectra (Fig.0\). 
To determine the goodness-of-fit of the phase-shifted models, we calculated as a function of the phase 
(j) and a scaling amplitude A, taking into account bin-to-bin correlations using the inverse Fisher matrix 
(Fig.E^). For this exercise we used the Al « 75 binning of the CBI power spectrum. The best-fit phase 
is 24° lb 33° with amplitude 0.94 relative to the fiducial model. The fiducial model is well within the 1 a 
(68%) confidence region (the difference in between the fiducial model and the best-fit model is 0.64 
for 2 df). The actual data and the best fit model are shown in Fig.0i;. This test shows that our data are 
entirely consistent with the model predictions, and that we can rule out (at w 3 o") a pathological model 
in which the EE oscillations are in phase with TT rather than out of phase. 

An alternate, and more physically motivated, way to look at the phase of the peaks in EE is to use 
fits to the fiducial model spectrum of the form 

/(/ + l)Cf ^ = (AJAso) {f{W/9o) + g{l9/eo) smikW/Oo)) . (2) 

This parameterizes the models in terms of two of the cosmological parameters discussed earlier, Ag 
and 6. The values of these parameters in the fiducial model are and 9q {Oq = 1.046). Changing 
6 scales the whole function, including the envelope, rather than just the phase. We now examine the 
variation of as these two parameters are changed, the other four cosmological parameters being fixed 
at their fiducial values. There is a minimum of near the fiducial model, with 9/9o = 1.02 it 0.04 and 
Ag/Aso = 0.93. (A second minimum in which the third polarization peak is shifted and scaled to fit the 
second fiducial peak is incompatible with the TT data.) This test also shows that the EE data strongly 
prefer the fiducial model, and demonstrates that the EE data alone have the power to place constraints 
on cosmological parameters. 

Tests for systematics. We have carried out a number of data quality tests to look for possible systematic 
contamination by foreground emission, residual ground emission, or other instrumental effects. We have 
found no evidence of significant residual instrumental or foreground effects after correcting for the point 
sources and projecting out the common ground spillover mode. 

Foreground emission is likely to have a different spectrum from the CMB, and ground contamination 
is frequency-dependent because it depends strongly on the baseline length in wavelengths, and thus 
shows up most on the shortest baselines at the lowest frequency. To look for these effects we estimated 
power spectra separately from the data taken in the lower and upper halves of our frequency band, i.e., 
26 to 31 and 31 to 36 GHz (fig. fTTl. We have compared the two spectra by using (including the 
bin-to-bin correlations). For EE and BB the measurements are dominated by thermal noise (rather than 
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sample variance) so the \^ results are valid. The results are: = 8.43 (7 df) for EE, and = 8.30 
(7 df) for BB. The probability of obtaining a larger by chance is 0.30 for EE and 0.31 for BB. 
The power spectra thus show no indication of strong contamination by foreground emission or residual 
ground emission. Note that for TT and TE the maximum likelihood error estimates include the effects of 
sample variance, and, because sample variance is correlated between the two frequency bands, a simple 
test is not valid. 

In addition to dividing the data into two frequency bands, we carried out jackknife tests in which we 
compared the following subsets of the data, (i) We compared all subsets of three of the four fields. This 
would indicate whether any of the fields is anomalous and is a good test for foreground contamination. 
No significant differences were found, and in particular the 14*^ field (which lies in the North Polar 
Spur region) was not anomalous, (ii) We compared all subsets of 12 of the 13 antennae. This would 
show up problems associated with particular antennae or receivers, (iii) We compared the TT spectra 
derived from the R and L antennae separately, to check for calibration discrepancies, (iv) We compared 
spectra estimated from the first and second halves of the data set, to check for effects based on season, 
distance from the primary calibrator, and other time-dependent parameters. None of these tests showed 
any significant differences between the data subsets. 

The DASI results increase our confidence that diffuse synchrotron emission is not a significant con- 
taminant in our EE spectrum. The fields that we have observed appear to be comparable to the DASI 
fields (Fig. Illl, and the DASI 95% confidence upper bound of 0.91 /xK^ on EE contamination should 
also apply to the CBI observations, which were made at higher / where the contribution of synchrotron 
emission is expected to be lower. 

Conclusions 

Our EE results are shown in comparison with the recent results from DASI and CAPMAP in Fig.[Sl We 
have detected the polarized CMB {EE) emission with high confidence (8.9 cr when foreground sources 
are ignored and 7.0 a when potentially contaminating sources are projected out), and we have also mea- 
sured the phase of the EE spectrum and shown that it is consistent with a phase-shift of vr relative to the 
TT, as expected if acoustic waves are the origin of the features in the TT and EE spectra on the scales 
of clusters of galaxies. The results from the CBI and DASI experiments are a powerful confirmation of 
the predictions of the standard model. The CBI continues to observe the polarized CMB emission, and 
we expect by the end of 2005 to have more than doubled the data set, leading to a decrease of over a 
factor two in the uncertainties of Q. 
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Supporting Material 

Modifications for polarization measurement The polarization observations reported in this paper 
were made between September 2002 and May 2004. Earlier tests of the technique were made using 
a single cross-polarized antenna The CBI was upgraded in 2002 to enhance its polarization 

capability. This involved replacement of the existing polarizers with new broadband achromatic po- 
larizers, replacement of the high-electron-mobility-transistor (HEMT) amplifiers with new lower-noise 
amplifiers, and reconfiguration of the antennae into a more compact array. 

The circular polarization mode (R or L) received by each antenna can be selected by changing the 
orientation of a quarter-wave plate in front of the low noise amplifier. The original CBI quarter-wave 
plates were replaced by achromatic DASI-style polarizers i67\\36[i that could be rotated under computer 
control so that the polarization in any antenna can be changed in < 5 s. An important design goal for 
polarization observations is to limit the polarization impurity: if an antenna does not receive pure R or 
L polarization the linear polarization measurements will be corrupted by an admixture of total intensity 
I. The fraction of / that appears in the RL or LR visibility measurement is called the leakage. This is 
a complex number (amplitude and phase) that must be measured for each antenna. The new polarizers 
reduced the leakages from ~ 5-15% to ~ 1-3%. The leakages are stable, and typically exhibit changes 
of < 0.2% over periods of a few months. 

In the first two years of operation of the CBI (2000 and 2001) we used sparse configurations of 
the antennae in order to cover a wide range of multipoles (300 < / < 3500). For the polarization 
observations we decided to concentrate on the multipole range 300-2000 in order to provide maximum 
sensitivity in the region where the CMB polarized EE signal is expected to be the greatest. For this 
reason we adopted the close-packed configuration shown in Fig.^This configuration provides the highest 
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concentration of short baselines possible with the CBI and provides an excellent match to the / range of 
the expected maximum EE signal. 

Data calibration and editing. During the observations and initial data analysis, we inspected each 
night's observations to look for instrumental and other problems. For the final analysis, we used auto- 
matic procedures to remove data with known problems (warm or unstable receivers, for example) and 
with higher than normal noise levels. This last check eliminated « 1% of the data that had been corrupted 
by clouds or instrumental problems. 

Amplitude calibration. For the CBI, the amplitude and phase calibration of the co-polar visibility 
data (RR or LL) was carried out using the same procedures as for the 2000-2001 observations (l?7l l. The 
refinement of the CBI flux density scale using the WMAP observations of Jupiter has been described 
in h24j . The uncertainty in the revised scale is 1.3% in flux density (2.6% in the power spectrum, Q). 
On most nights one or more of the primary calibration sources Tau A (the Crab Nebula), Jupiter, Saturn, 
and 3C 274 was observed. All of these gave consistent results, except for 3C 274: we found that the 
flux density of 3C 274 declined by 7% over the period of these observations. The model for 3C 274 was 
adjusted to take this secular variation into account before the final calibration of the data. 

The majority of CBI data are calibrated using measurements of Tau A and Jupiter, as in (24). When 
none of the primary calibration sources was available we used secondary calibration sources, such as the 
variable quasar J1924— 293, for which we obtained flux densities by interpolating from adjacent days 
calibrated against the primary calibrators. 

Polarization calibration. The calibration measurements on the co-polar baselines yield complex gain 
factors for each antenna. These gain factors are sufficient to calibrate the cross-polarized baselines {LR 
and RL) except for an unknown phase difference between the R antennae and the L antennae, equivalent 
to an unknown rotation of the plane of linear polarization We determine the unknown L — R 

phase difference by observations of a strong, polarized calibration source, Tau A, for which we assume 
the polarization position angle (i?- vector) is —27.6°. This value was derived L65..66) by comparison of 
CBI and Very Large Array iSTh observations of 3C 273 and 3C 279, both of which vary but are observed 
regularly with the CBI at 26-36 GHz and the VLA at frequencies straddling the CBI band (22 GHz and 
44 GHz) {68). It is close to the angle measured with other instruments at lower frequencies. The L — R 
phase difference is very stable unless receivers or cables are modified, so on nights when no measurement 
of Tau A was available we used the average of all the Tau A measurements. 

Leakage measurement We measured the instrumental polarization leakage factors on each night when 
either Tau A or Jupiter could be observed. These observations were made at a number of different 
parallactic angles, by rotating the CBI platform, to enable the source and instrumental polarization to be 
separated. The instrumental leakage was found to be in the range l%-3% on most baseline-channels, 
with a few baseline-channels showing leakages as high as 5%. We determined that leakage did not vary 
significantly across the field of view by making observations of Tau A at a number of offset positions. 
This low level of instrumental polarization and our strategy of rotating the deck so that many antenna 
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pairs contribute to the same (u, v) point ensures that any instrumental polarization in our final data set is 
negligible (< 1%). We have therefore ignored the instrumental polarization in the present analysis. 



Noise calculation. In order to ensure that the ground contamination was identical in each of a set 
of six pointings, we deleted all visibility samples that did not have counterparts observed at the same 
hour angle (within the tolerance of the integration time, 4.2 s) in all of the six fields. After selecting 
matched data points in this way, we calculated the noise from the scatter of the visibility measurements. 
As an error in the noise estimate will bias the final power spectrum estimate, it is important to obtain 
an accurate estimate of the noise in the data. In one scan, comprising observations of 6 fields, we 
record m = 1, . . . , M (M varies, but is usually about 37) data points (complex visibilities) for each of 
n = 1, . . . , fields (A^ = 6) . The observed visibility Vnm is related to the true visibility Xn and the 
ground contribution g-m by 

^nm — Xn ~\~ Qm ~l~ 1"nmi (3) 

where r„m is the noise in the measurement. Note that Xn is the same for all m (we do not change the 
baseline length or orientation relative to the sky during the scan), and gm is the same for all n (the ground 
contribution is assumed to be the same in each field, i.e., constant for the duration of the scan at a given 
elevation and azimuth). An estimator for the noise variance is 

N M 

(M- 1)(A^- 1) ^ "™ 

where is derived from Vnm by subtracting the mean of the N measurements from all N obtained 
at each time m, and the mean of the M measurements from all M obtained on each field n. Our best 
estimator is the average of the two estimates obtained by treating the real and imaginary parts of the 
visibility separately. We obtained a single noise estimate for each baseline-channel that applies to a 
whole 18 min scan. The variance of the estimator is 

^■"■f^' = (A.-lKiV-l) - 

The uncertainty in the noise estimate in each scan is small enough that noise bias is not a concern 
in the present observations. Scans with rms noise more than three times that expected for normal system 
temperatures were deleted; in most cases the high noise was due to clouds. 



Power spectrum estimation The principles of estimating polarization power spectra from interfer- 
ometer visibility measurements are described by Wh . To process the CBI data, we have extended the 
gridding-based procedure used in our earlier work (l49l l to deal with mosaicked polarization observations. 
A given correlator output sample, or visibihty, can be one of the four polarization products RR, RL, LR, 
or LL. These can be related to the fundamental CMB polarization modes T (temperature), E, and B 
(polarization) through the expressions given in Equations 3 and 4 of {37). The resulting power spectra 
are decomposed into the six possible covariances TT, EE, BB, TE, TB, and EB. Note that because 
the CBI measures circular polarization products, which are orientation independent (depending only 
on the handedness of the wave polarization), the CBI (or any interferometer using circularly polarized 
receptors) is sensitive to the E and B modes directly. This simplifies the power spectrum analysis. 
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The co-polar RR and LL visibilities are gridded together into an effective RR estimator ({LL) 
and {RR) are identical in the absence of circular polarization) as in (1?^ . while the cross-polar RL 
and LR visibilities are gridded together, after conjugating and reflecting the LR visibilities in the uv- 
plane, into cross-polar estimators A/j^ using the same gridding kernel as the co-polar data. The covari- 
ance matrix elements are computed for the cross-polar estimators using a modified operator Pjil{v) = 
F{v) e*^(^~'^) where P is defined in Equation 12 of ( 1?^ . ip is the on-sky parallactic angle of the CBI 
receivers (Equation 2 of {36 )) and x is the wave- vector angle corresponding to the uv point v (Equation 3 
of L37)). The band powers derived from the likelihood analysis are then {q%,B = 1, . . . , N^}; the dif- 
ferent covariance products S = TT, EE, BB, TE, TB, EB can have different numbers and locations 
of bands. Point sources are handled in the same manner as in {49), with the option of projecting out the 
sources from the RL and LR parts of the covariance. The new scanning procedure required the addition 
of a scan projection matrix C^'^^'^ constructed by building a "noiselike" matrix as in Equations 32 and 35 
of with the covariance elements E^ki = 1 if visibilities k and k' are from the same scan (otherwise 
zero); this is then projected out with a pre-factor gscan in the same way as the point sources are. Details 
of this modified procedure will be given in ( IS?b . 

The maximum likehhood estimation of the spectrum from the gridded estimators is done on the CITA 
McKenzie cluster dS^t which consists of 256 nodes with two 2.4 GHz Intel Xeon processors and 1 GB 
of memory per node. The matrix operations are done using the SCALAPACK library STUii . From an 
initial guess of the spectrum, we iterate to the maximum likelihood solution using the Newton-Raphson 
method. One modification to the procedure used in i53l provides a significant improvement. Rather than 
use the standard approximation to the curvature, with which the number of expensive matrix operations 
is proportional to the number of bins in the spectrum, we use an approximate curvature that requires 
only a single matrix inversion dZTT i. Using 32 nodes per mosaic, with 10'^ estimators per mosaic, this 
decreases the time per iteration from about an hour to one minute, without changing the solutions. The 
total time for the spectrum to converge, once the estimators and correlations are read into memory, is 
about 10 minutes, and is virtually independent of the number of bins in the spectrum. 
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Table 1: CBl Band Powers in fjK^. The broad first band is not sensitive over the 0-600 range, but with 
finer binning we have a detection in a band near 400. 
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Table 2: Cosmological constraints from the WMAPl only, CBI pol plus WMAPl, and CBI pol plus 
CBIext plus WMAPl data compilations for an assumed Q,tot = 1-0. Relaxation of this constraint opens 
up the tight uncertainties on Hq and O^. We included weak external priors on the Hubble parameter 
(45kms~^ Mpc~^ < Hq < 90 kms~^Mpc~^) and the age of the universe (to > lOGyear). The 
flatness prior has the strongest effect on the parameters by breaking the geometrical degeneracy and 
allowing us to derive tight constraints on Hq and 0,^,- The top six parameters are those used in the 
Markov chain evaluations and the distributions of the bottom six are derived from the same chains. The 
uncertainties are given as 68% confidence intervals. 
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Figure 1: The Cosmic Background Imager. For the polarization observations, the 13 90-cm Cassegrain 
antennae were arranged in this hexagonal close-packed configuration on the rotating, alt-az mounted plat- 
form, with six adjusted to be sensitive to right-hand circularly polarized radiation and seven to left-hand 
circularly polarized radiation. By correlating the signals from the antennae in pairs, 78 interferometer 
baselines are obtained ranging in length from 1.0 to 3.5 m. 
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Figure 2: The sky coverage of the four fields imaged by the CBI in polarization. The grey-scale shows 
the noise-level achieved in total intensity, /, in the observations reported here. Three of the fields have 
been mapped with 36 separate pointings, whereas the fourth, 20^^, has been mapped more deeply but 
in only six pointings. The modulation of the sensitivity by the CBI primary beam is apparent. The ap- 
proximate centers of the four fields are: 02'^49'^30^ -02°52'30"; 08^47'^30^ -02°47'30"; 14^45'"30^ 
-04°07'30"; 20'^49'^30^ -03°30'00" (J2000 right ascension and dechnation). 
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Figure 3: Location of CBI and DASI fields in relation to the Galaxy. The sky image is the Ka-band 
synchrotron map derived from WMAP first year data (.72.) . Galactic longitude increases to the left, with 
zero in the center of the image. 
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Figure 4: Images of the 14^ field mapped by tiie CBI in Stokes parameters /, Q, and U (Stokes V, circular 
polarization, is not measured and is assumed to be zero). Color is used to represent intensity, with the 
same scale in each Stokes parameter. In these images the contaminating effects of ground radiation 
and foreground emission have not been removed. The total intensity, /, image {left) is dominated by 
CMB emission (modulated by the instrumental point-spread function); some foreground point sources 
are visible {red spots). The linear polarization, Q and U, images {center and right) are dominated by 
instrumental noise and ground pickup. Ground pickup, which with our observing strategy should be the 
same in each pointing at the same declination, gives rise to a pattern that repeats at intervals of 3 min in 
right ascension. 
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Figure 5: The effect of lead minus trail differencing. Here the data presented in Fig. |4] have been dif- 
ferenced: each visibility measurement has had the corresponding measurement on a field 9 min later in 
right ascension subtracted. Because the ground pickup is very similar for both measurements, ground 
emission cancels out in the difference. In the resulting images foreground point sources may appear pos- 
itive or negative in I depending on their right ascension. The Q and U images show that ground pickup 
has been removed with high accuracy. 
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Figure 6: Power spectra of CMB polarization from the CBI measurements. The four panels show total 
intensity power spectrum TT, grad polarization mode power spectrum EE, curl polarization mode power 
spectrum BB, and cross-spectrum TE. Numerical values are given in Tabled The black curve is the 
theoretical ACDM model using a power law for the primordial spectral index which best fits the WMAP, 
CBI, and ACBAR CMB data (57). The predictions of this model for the CBI bands using the CBI 
window functions are indicated by the stars. 
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Figure 7: Sensitivity of the CBI EE power spectrum to the phase of the acoustic oscillations. (A) The 
EE spectrum predicted by the fiducial model (black line) with a variety of "phase shifted" spectra with 
similar envelopes calculated as described in the text. (B) Goodness-of-fit (x^) for the model as a function 
of two parameters: the phase shift (horizontal axis) and an overall scaling (vertical axis); the point at (0,1) 
indicates the fiducial model corresponding to a standard ACDM model. Contours are at 1, 2, and 3 a 
intervals (i.e., Ax^ = 2.30,6.17, and 11.8 for 2 df). (C) Comparison of the fit of the fiducial model 
(black line) and the minimum-x^ phase-shifted model (red line) with the CBI data points (blue) and the 
band powers predicted by the models (black and red stars). 
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Figure 8: Comparison of EE measurements from CBI, DASI OH, and CAPMAP OB- The fiducial 
model curve is the same as in Fig.[6l The asterisks show the predictions of the fiducial model for the CBI 
bands. 
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Figure 9: Normalized band-to-band correlations for the 7 bands shown in Table[l] 
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Figure 10: Band-power window functions for the CBI polarization observations with the bands shown in 
Table[2 The expected value of band power for a given model spectrum, C/, is J2i [^(0/^] ^(^ + l)C'i/27r. 
The total area under each window function is equal to unity. The band powers for each spectrum also 
contain contributions from the the others (EE band powers are affected by changes to the model TT and 
TE as well as EE, for example), but for the CBI these cross-polarized window functions are small, with 
peak amplitudes of no more than a few percent of the same-polarized window functions. 
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Figure 1 1 : CBI polarization power spectra obtained from low- and high-frequency channels. Red points: 
26-31 GHz; blue points: 31-36 GHz. The points have been offset horizontally in / for clarity. The four 
panels show total intensity power spectrum TT, curl-free polarization mode power spectrum EE, curl 
polarization mode power spectrum BB, and cross-spectrum TE. The fiducial model curve is the same 
as in Fig. 6. 
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